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5:^ ' Abstract 

We propose discrete mappings of second order that have a discrete analogue of 
Lyapunov function. The mappings are extensions of the integrable Quispel-Roberts- 
Thompson (QRT) mapping, and a discrete Lyapunov function of the mappings 
is identical to an explicit conserved quantity of the QRT mapping. Moreover we 
■ can obtain a differential and an ultradiscrete limit of the mappings preserving the 

, existence of Lyapunov function. We also give applications of a mapping with an 

Q \ adjusted parameter, a probabilistic mapping and coupled mappings. 
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1 Introduction 



Among discrete mappings, integrable ones have a rich mathematical structure 
such as explicit exact solutions, conserved quantities, symmetries, and so on. 
For example, most discrete soliton equations belong to the discrete KP hierar- 
chy, have explicit N soliton solutions and have an infinite number of explicit 
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conserved quantities[l,2,3]. The discrete Painleve equations have a direct re- 
lation to the affine Weyl group and their mathematical structure of explicit 
solutions gives a hot theme[4,5]. Integrable systems reflect some aspects of 
essential mechanism of real nonlinear phenomena in conserved systems and 
studies on the systems are important subjects. 

On the other hand, general discrete systems are complex and do not exhibit 
such explicit mathematical structure as in the case of integrable systems. How- 
ever they reflect wider nonlinear phenomena and show a rich structure of so- 
lutions, which can not be easily grasped by the integrable system theory. Thus 
extending the theory of integrable systems is an important theme. 

In this paper, we show a link between integrable and nonintegrable mappings. 
In section 2, we propose a certain type of discrete mappings that possess a dis- 
crete analogue of the Lyapunov function, which increases or decreases mono- 
tonically and converges to a certain value in the course of time evolution[6]. 
We call this function a discrete Lyapunov function for short. The mappings are 
obtained by extending a class of second-order discrete mappings, the Quispel- 
Roberts- Thompson (QRT) mappings[7,8,9,10]. The QRT mappings are inte- 
grable and have an explicit conserved quantity that remains constant in the 
course of time evolution. The conserved quantity of QRT mappings becomes a 
discrete Lyapunov function in the corresponding extended mappings we pro- 
pose and a trajectory of solution gets close to an attractor which is deflned by 
the discrete Lyapunov function. 

In section 3, we give a differential and an ultradiscrete limit of the mappings[ll,12,13]. 
The ultradiscrete equation is deflned by using ± and max (or min) operators. 
The ultradiscretizing procedure was discovered in the integrable system the- 
ory and various discrete soliton equations can be ultradiscretized through a 
certain type of transformation of variables and parameters. The differential 
equations and the ultradiscrete equations obtained from the mappings have 
also an explicit Lyapunov or discrete Lyapunov function derived from that of 
the mappings. 

In section 4, we give some applications of the mappings. Preserving mono- 
tonicity of the discrete Lyapunov functions, we can construct three types of 
variant of the mappings, a mapping with an adjusted parameter, a probabilis- 
tic mapping and coupled mappings. In section 5, we give concluding remarks. 

Before going into the main subject, we make a short introduction to the QRT 
mapping and to the limiting procedures. The symmetric version of QRT map- 
ping is deflned by the following second-order discrete equation[7,8], 
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where fi{x) ~ fsix) are defined by 
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with symmetric 3x3 matrices A and B, 
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Note that 'x' is an outer product. If we define hn — h{xn-i,Xn) by 

N{x,y) 



h{x,y) 



D{x,yy 
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where 



N{x, y) = [x^ xl]A 



D{x,y) = [x^ xl]B 



(5) 



we can easily show = hn- Thus hn is a conserved quantity of (1) and the 
QRT mapping is integrable. 

We can obtain a differential equation from the QRT mapping[14]. For example, 
let us consider a QRT mapping 



Xn+l 



Xji + + '■'2 

{CsXn + l)Xn-l 



(6) 



with a conserved quantity 
1 



h(^Xn—liXn) — 



Xn-l+Xn+Cl(Xn-l+Xn)+C2+C3Xn-lXn{Xn-l+Xn)). (7) 



Xn—l-^n 

If we use a transformation including a new parameter S defined by 



Xn 6 



y{n5) 



Cj — S Cj, 

a series expansion of (6) at S ^ gives 

/ = c^e-y + C2e-'^ - cse^ + 0{6). 



(9) 
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Thus we obtain a differential equation 

y" = c^e-y + C2e-2^ - cge^ (10) 

as a limit equation of (6). This equation is also integrable since it has a con- 
served quantity derived from /i, 

h{y,y') = liml(/i(e^(*-'^),e^W) - 2) = {y'f + 20^6-^ + c^e-^^ + 2c3e^. (11) 

We can also obtain an ultradiscrete mapping from the QRT mapping[ll,12,13]. 
Consider (6) again and use a transformation of variable and parameter includ- 
ing a new parameter e defined by 

Xn = e^"/^ c,- = e^^''. (12) 

Then we obtain an equation on X„, 

^Xn/e _j_ g(X„+Ci)/c _|_ qC'z/s 
= + £ log ^(x„+C3)/. + l ■ ^^^^ 

Taking a limit £ ^ -|-0 and using a formula 

^lim^ £ log(e^/" + e^'^ + ■•■)= max(A, • ■ • ) , (14) 

we obtain 

Xn+i = -Xn-i + max(2X„, X„ + Ci, C2) - max(X„ + C3, 0). (15) 
This equation is also integrable since it has a conserved quantity derived from 

K 

H{Xn-i, Xn) = ^hm £ log /i(e^"-l/^ e^"/^) 

= max(X„_i + X„ + max(X„_i, X„) + C3, 2X„_i, 2X„, (16) 
max(X„_i, Xn) + Ci, C2) - Xn-i - X^. 



2 Derivation of mapping 



In this paper, we discuss discrete mappings which have the following form, 

^ fl{p^n) •^n—lf^ip^n) Oin-^n— l(^n ^oo) fl7) 

where /i(a;) ~ fsix) and hn are the same as those defined in (2)~(5) and h^o 
is an arbitrary constant. Let us consider a solution x„ for n > from initial 
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data xo and xi without loss of generality. Condition for q;„ is shown below. 
Noting 

hn+l ~ hn 

(f2(Xn) - Xn-lfsiXn)) - ifl(Xn) " Xn-lf2(Xn)) 



D{Xn-l, Xn)D{Xn, Xn+l) 



and using (17), we obtain 



(18) 
(19) 

D{Xn-l, Xn)D(Xn, Xn+l) 

Therefore, if e is a small positive constant and < «;„ < 1 — e holds for any 
n, then hn decreases or increases monotonically and converges to hoo- This 
condition is rewritten by 



where 



(Xn+l — Xn-l)"^ 



1- ar, 



e < a. 



(^n+l ~ Xn-x)'' 



D{Xn-l, Xn)D{Xn, Xn+l) 



< 1. 



(21) 



If this is satisfied, any trajectory of solution gets close to an attractor /i„ = h^^, 
in a phase plane and hn becomes a discrete analogue to Lyapunov function of 
a differential equation. 

If we solve an inequality (21) on q;„ with (17), a general condition of is 
expressed by Xn-i and ,t„. However, since (17) itself includes and a sign of 
D{xn^i, Xn)D{xn, Xn+l) depends on a:„_i, a:„ and x^+i, the condition becomes 
quite complex. Therefore, we choose special matrices A or in (3) and propose 
a sufficient condition of an to make the mapping (17) rational. 



If we choose a matrix B as 



B 





1 



(22) 



then fj, D and h become as follows: 

fi{x) = aoix"^ + aux + ai2, 
f2{x) = -{aoox'^ + aoix + 002), 
Mx) = 0, 
D{x,y) = l, 

h{x, y) = aoox^y'^ + aoixy{x + y) 

+ 002(2;^ + y^) + aiixy + 012(2; + I/) + a22- 



(23) 
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In this case, (1) is of class I after Refs. [9,10] and (17) becomes 



h[Xn) - OLn[hn - /ioo) 

The condition (21) is reduced to 

f fljXn) - 2Xn-lf2{Xn) ^ ^ ^ , . 

The latter inequahty gives 

Oin{{h{Xn) - 2x„_i/2(x„))^ + 2f2{Xn){K " ^oo)) < f2{Xnf + ^liK " ^oo)^- 

(26) 

Then we can derive a sufficient condition for q;„, 

' ^ - ifliXn) - 2Xn-lf2ixl)r + ^(Xn)^ + (h^ - h^f ^^^^ 

using 2 /2(a;„) (/;,„- /loo) < f2{xnf + {hn-hoof and < al{h„-h^y. Figure 1 
shows solutions to (24) plotted in the phase plane using q;„ defined by 

^ f2{XnY , . 

'''' ~ ^ ifM - 2x^_Mxr,)r + Mxr^r + (h^ - h^r ^ ^ 

with an additional parameter p satisfying < p < 1. Note that f2{x)'^ = 
(x^ + 2x + 3)^ > 4 and the condition e < a„ in (27) is satisfied in this case. 
The only fixed point Xn ~ —0.1649 exists and does not depend on p. Figure 1 
(a) shows a solution to (24) with p — 1 and {xo,xi) = (3,2). The trajectory 
rapidly gets close to an attractor hn — h^o in the outer region. Parameters of 
Fig. 1 (b) arc all the same as of Fig. 1 (a) other than (xo,Xi) = (—0.2, —0.3). 
This trajectory moves in the inner region of the attractor and also approaches 
it rapidly 

If /t„ = in (19), hn+i = hoo holds and points (a;„_i,a;„) for 2 < n are all 
on the attractor hn = h^o- If /^n = 1, hn+i = hn holds and all points are on 
an intcgrable trajectory defined by initial data and do not get close to the 
attractor. Therefore, a closing rate of solution to the attractor depends on 
that is, «„. From the relation (20) between and we see that the closing 
rate is smaller if q;„ is smaller under (27). Figure 2 shows solutions in the case 
of p = 0.02 with the same parameters as in Fig. 1 other than p. Since is 
smaller than that of Fig. 1, the trajectory gets close to the attractor more 
slowly. Figures 3 (a) and (b) show an evolution of hn for Figs. 2 (a) and (b) 
respectively. In both figures /i„ converges to hoo (=5) monotonically. 
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-4 -2 2 4 




(a) 



(b) 



Fig. 1. Solutions to (24) with q:,„ defined by (28) for 
(aoo,aoi,ao2,aii,ai2,a22) = (1,2,3,1,1,0), hoo = 5 and p = 1. Points (x„_i,a;„) 
(1 < n < 1000) in the phase plane are plotted, (a) (xo,xi) = (3,2), (b) 
{xo,xi) = (-0.2,-0.3). 




-4 -2 2 4 



(a) 



(b) 



Fig. 2. Solutions to (24) with defined by (28) for p = 0.02. Parameters other 
than p are the same as in Fig. 1. 
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Fig. 3. Evolution of hn for Figs. 2 (a) and (b). 
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3 Differential and ultradiscrete equations obtained from discrete 
mappings 



In this section, we consider a special case of the mapping (17) by choosing B 
in (3) as 





B = 



1 




(29) 



Then fj, D and h are given by 

fi{x) = -x{ao2x'^ + ai2X + 022), 

f2{x) = 0, 

/3(a;) = x{aQQX^ + a^^x + 002), 
D{x,y) = xy, 

H^: y) = — {aoox^y'^ + aoixy{x + y) 
xy 

+ 002(2:;^ + y'^) + a^ixy + 012(2: + y) + 022). 
In this case, (1) is of class II after Refs. [9,10] and (17) becomes 

gi{Xn) h^gi{Xn) + Kxl_^g^{Xn) + 



(30) 



Xn+l 



where 



Xn-ig^{Xn) Kgi{Xn) + h^xl_-^g2,{Xn) + /?n ' 

gi{x) = -fi{x)/x = aQ2X^ + ai2X + 022, 
^3(2;) = f^{x)/x = aoox^ + aoix + 002, 

Xn—lXn9l iXfi)g^ (-^ri) 



(31) 



Otr, 



Note that the mapping 



hoc{gi{xn) +xl_igs{xn)) +/3n 

_ gi{xn) 

Xn+1 



Xn-igsiXn)' 

is the original QRT mapping of class II. Assume 



(32) 



(33) 



(i) ttij > for any i and j, 

(ii) if X > and y > 0, then gi{x) > 0, 513(2;) > and h{x, y) > 0, (34) 

(iii) /loo > and /3n > 0. 

Then, if Xn-i > and x„ > 0, we obtain Xn+i > 0. Therefore, x„ from positive 
initial data is always positive. If we assume that x„ is always positive and that 
(34) is satisfied, we can prove (21) is also satisfied. The proof is complicated 
and is omitted here. Under this assumption, any trajectory of solution other 
than a fixed point gets close to the attractor hn = /ioo in the phase plane for 
any positive Figure 4 shows examples of solution to (31). 



8 



1 2 3 4 5 6 

Fig. 4. Solutions to (31) from {xq,xi) = (5,1) and (2,1) with 
(aoo,aoi,ao2, an, 012,022) = (1,0,0,0,1,2), ^00 = 6, = 100. Points (xn-i,x„) 
(1 < n < 1000) in the phase plane are plotted. 

Moreover, there is a remarkable feature about the mapping (31). Differential 
and ultradiscrete (piecewise-linear) equations are obtained by two kinds of 
limiting procedure. 



3.1 Differential equation obtained from (31) 



With regard to a differential equation, we use the following transformation 
including a small parameter 5. 





— gJ/(n5+Co) 






O02 


= C02, 


On 


= -2co2 + S^Cii, 




= ^ Coo, 




= S'^ Coi, ai2 




= (5^ ^00, 


Pn 


^6/ p. 



(35) 



Substituting (35) to (31) and taking a limit 5 — > 0, we obtain 

/ = — f-Cooe'^-Coie^ + Ci2e-f + C22e-2^) -2co2Pe''^y'(h{y,y')-h^), (36) 
C02 V / 

where h{y, y') is defined by 

h{y, y') = Co2(t/')' + cooe'^ + 2coie^ + Cn + 20126"^ + 0226-^^. (37) 

Note that /i(7/, 7/') can also be obtained from the discrete Lyapunov function 
hn — hi^Xn—i, Xyi) as 

h(y, y') = hm ^/i(e^«"-^)''+«°), e^'^^^+^o)). (38) 
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(a) (b) 

Fig. 5. Solutions to (36) for (cqo, cqi, C02, cn, C12, C22) = (2,1,1,0,1,1), hoo = 8, 
(3 = 0.2 and (a) y(0) = 1, y'(0) = 0, (b) y(0) = y'(0) = 0. Both solutions are 
calculated numerically by the fourth-order Runge-Kutta method. 

Using (36), we can show 

jfi - hoo) = -4clJe'y{yf{h - h^). (39) 

This means h converges to and h is a Lyapunov function to (36). Figure 5 
shows examples of solutions in a phase plane {y,y'). The trajectories become 
a spiral shape approaching to a limit cycle defined by h — hoo- 



3.2 Ultradiscrete equation obtained from (31) 



We use the following transformation including a new parameter e to ultradis- 
cretize (31), 

a,,• = e^-/^ a:„ = e^"/^ /i^o = e^°°/^ /3„ = e^"/^ (40) 

Substituting these relations into (31) and taking a limit e +0, we obtain 
an ultradiscrete equation, 

Xn+l — Gi{Xn) — Xn^i — Gz{Xn) 

+ max(//^ + Gi{Xn), + 2X„_i + G^{Xn), Br,) (41) 
- max(//„ + ^(X^), Hoo + 2X„_i + G^{X,,), S„), 

where 

G'i(X) = max(Ao2 + 2X, A^^ + X, A22), 
G^{X) = max(Aoo + 2X, Aqi + X, A02), 

Hn = H{Xn-l,Xn) , , 

H{X,Y) =m£ix{Aoo + 2X + 2Y, Aoi+X + Y + msix{X,Y), ^ ' 

A02 + 2 max(X, Y),Au+X + Y, A^2 + max(X, F), ^22) 
-X-Y. 
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Fig. 6. Solutions to (41) with (^00,^01,^02,^11,^12,^22) = (13,71,11,17,43,23), 
i?oo = 100. (a) ixo,xi) = (800,700), (b) (-13,-12). 



Note that = corresponds to Aij — > — 00 and the terms including Aij in 
(42) are eliminated from max functions. Figure 6 shows examples of solution 
to (41). Trajectory is exactly on the attractor Hn = H^o after a finite number 
of time steps in both figures. Figure 7 shows an evolution of points in a region 
of the phase plane by the mapping with the same parameters in Fig. 6. Points 
{x,y) — (5i, 5j) (—40 <i,j< 40) are taken as initial data. All points in this 
region are mapped to the attractor after 7 steps. Note that the only fixed 
point is (—14, —14) and it is not shown in this figure. 

The polygon in the phase space defined by H^, = H^o is always an attractor 
of (41). However, additional attractors of integrable trajectory often appear 
depending on the parameters of (41). To explain this phenomenon, let us 
consider a relation a = 1 + 6. This relation gives an inequality a > b. If 
we use a transformation a = c^/^, b = e^^^, a relation between A and B is 
A = £log(l + e'^/^) and A > B also holds. However, A = max(0, B) is obtained 
after the limit £ — > +0 and A> B holds, especially A — B ii B > 0. This type 
of change of inequality often occurs in the ultradiscretization and also for the 
convergence condition (21). It is a reason why additional attractors appear for 
(41). 



4 Some applications 



We give some applications of mappings we proposed in the previous sections. 
Mappings in the applications also have an explicit discrete Lyapunov function. 
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Fig. 7. Evolution of points = (5z,5j) (—40 <i,j< 40) by the mapping (41). 

All parameters are the same as in Fig. 6. 



4-1 Parameter adjustment 



If a solution satisfies (21), /i„ converges to h^o and the convergence speed 
depends on Therefore, we can easily control the local speed by adjusting 
«„. Consider the case that p in (28) depends on n. Even in such a case, if p 
satisfies < p < 1, /i„ converges to h^o monotonically. For example, let us 
assume two constants cq (0 < cq < 1), ci and define p by 

p = Co + (1 - Co) sin^(ci(/i„ - /loo))- (43) 

Figure 8 (a) shows an example of solution in a phase plane. The parameter p 
depends on hn — hoo and p is approximately equal to Co near the region where 
sin(ci(/i„ — /loo)) ~ 0. Therefore, if we use a small cq, the convergence speed 
becomes slow around the region and a trajectory of solution there becomes 
similar to an integrable one. Figure 8 (b) shows an evolution of hn of the 
solution in Fig. 8 (a). It changes stepwise and converges to h^o finally. 
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Fig. 8. Mapping (24) with of (28) adjusted byp = co + (l — co)sin^(ci(^„ — ^oo))- 
Parameters and initial data are the same as in Fig. 1 (a), cq = 0.00004 and c\ = 0.2. 
(a) Solution in a phase plane, (b) Evolution of 



4-2 Probabilistic mapping 



Wc can easily introduce a probabilistic parameter to preserving the conver- 
gence. This application is based on a similar idea to the parameter adjustment. 
Consider the case that p in (28) is a probabilistic parameter. Even in such a 
case, if p satisfies < p < 1, /i„ converges to hoo monotonically. As an extreme 
example, let us assume that a value of p is or positive c with a probability 
defined by 

P(p = 0) = l-^, P(p = c) = ^, (44) 

with a positive integer M. The mapping (24) becomes integrable at time steps 
of p = {an = 0) and hn becomes constant for those steps. However, hn 
converges to hoo as n — > oo since the ratio of time steps of p = c is nonzero 
finite. Figure 9 (a) shows a solution and Fig. 9 (b) shows an evolution of hn- 



4-3 Coupled mappings 



The solution x„ to (31) satisfying (34) is always positive when positive initial 
values are used. Moreover, the only condition required for Pn in (34) is /3„ > 0. 
Therefore we can easily couple some of mappings (31). The general form of a 
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-4 -2 2 4 5 1- ---- ^ , ^ 

Xn-l 500 1000 

(a) (b) 

Fig. 9. Probabilistic mapping (24) with a„ of (28) using probabilistic p defined by 
(44) with c = 0.5 and M = 30. Other parameters are the same as in Fig. 1 (a), (a) 
Solution in a phase plane, (b) Evolution of 

system of x^^^ , x^^^ , • • • , x^j^^ is defined by 



X. 



fa) gP(^g') 



X 



/ii^'^^?'\a;(/')) + h^ixliL^'g^^ix'i^) + (30){n, x^l„ • • ■ , x^, • • • , xi""^) ' 

(45) 

where 

h^^^^Hx^^l^x^^^). (46) 

Note that gf^^ gl^^\ h^^ and /S^-') can be defined independently for every j and 
(3^^ is an arbitrary positive function. 

Figure 10 shows a simple case, coupled two mappings (M = 2) where 

(1)^ A (2)/ \ (1)/ \ (2)/ N 

9i{x)=g\'[x), 9l{x)=gl'{x), 

h('\x,y) = h('\x,y), = (47) 

= 1000, = exp(100(/iW - /i^^^)). 

In this case, x^^ is not affected by x^\ but x^^ is by x^^^ through If 
initial data of satisfy /iq^^ > /i^^ , then /3^'^^ > and the mapping for x^^^ 
is similar to an integrable one initially. Therefore the convergence speed of h^^^ 
to is slow initially. However, x!^^ satisfies h^^^ ~ for n ^ 0, then P^'^^ 
becomes small and the convergence speed of h^^ becomes larger. Figure 10 
(a) and (b) show an evolution of h^^ and of h!^^ respectively. Figure 11 (a) 
and (b) show solutions x!^^ and x!^^ in a phase plane respectively. 
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Fig. 10. Evolution of the discrete Lyapunov function (a) hn^ and (b) /i^f ^ . Common 
parameters to define g[^\ hn^ are (000,001,002,011,012,022) = (1,0,0,0,1,2), 



"-00 



I loo 
(2) _ 



5, = 1000, = exp(100(/i^/^ - /i^^)), x, 



,(1) 



.(2) 
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Fig. 11. Solutions in a phase plane for (a) Xn^ and (b) Xn^. Mappings and initial 
data are the same as in Fig. 10. 



5 Summciry and concluding remarks 



Wc propose second-order discrete mappings as an extension of the QRT map- 
pings. They have an explicit discrete Lyapunov function and the asymptotic 
behavior of a solution is clearly understood. The function is the same as the 
conserved quantity of the corresponding QRT mapping. This fact suggests 
that we can extend an integrable system to nonintegrable one together with 
its characteristic structure. 

Moreover, a differential and an ultradiscrete equations corresponding to the 
mapping can be derived by taking a limit of a parameter. Both equations have 
also an explicit Lyapunov function. Thus we give a set of discrete, differential 
and ultradiscrete equations of which behavior of solutions obey the same struc- 
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ture. There are many examples of this kind of set among integrable equations. 
It suggests that there may be some mathematical structure in the equations 
proposed. Note that although it is not so difficult to construct a differential 
equation with an explicit Lyapunov function, it is difficult to discretize or 
ultradiscretize the equation together with its Lyapunov function. 

We show some applications of the discrete mappings. We give a mapping 
with an adjusted parameter, a probabilistic mapping and coupled mappings 
preserving the discrete Lyapunov function. As shown in the previous section, 
the existence of the function makes it possible to construct applicable systems 
of which behavior is easy to control. 

Lastly we list below some future problems. 

• It is clear that the mapping proposed is nonintegrable since they arc not 
reversible. However, more detailed study may reveal a relation between the 
mapping and its integrable correspondence. 

• There may be another characteristic quantity of the mappings. Utilizing the 
quantity together with the Lyapunov function, we may be able to give more 
information about solutions. 

• Though we concentrate on modifying the QRT mapping, a similar procedure 
may be applied to other integrable mappings. Generalizing the method used 
is an interesting problem. 

• We propose the applications of mappings only as demonstrative examples. 
We should find an application as a model of more realistic phenomenon. 
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